Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 14 September 2012 



(MN WI^ style file v2.2) 



The ballistic transport instability in Saturn's rings I: 
formalism and linear theory 



(N 
O 

m 



o 

in: 



> 

o; 



Henrik N. Latter^*, Gordon I. Ogilvie^f, Marie Chupeau^'^J 

1 DAMTP, University of Cambridge, CMS, Wilberforce Road, Cambridge CB3 OWA, UK 
^ ENS Cachan, 61 Avenue du President Wilson, 94230 Cachan, France 



ABSTRACT 

Planetary rings sustain a continual bombardment of hypervelocity meteoroids that 
erode the surfaces of ring particles on time scales of 10^-10^ years. The debris ejected 
from such impacts re-accretes on to the ring, though often at a slightly different orbital 
radius from the point of emission. This 'ballistic transport' leads to a rearrangement of 
the disk's mass and angular momentum, and gives rise to a linear instability that gen- 
erates structure on relatively large scales. It is likely that the 100-km undulations in 
Saturn's inner B-ring and the plateaus and 1000-km waves in Saturn's C-ring are con- 
nected to the nonlinear saturation of the instability. In this paper the physical problem 
is reformulated so as to apply to a local patch of disk (the shearing sheet). This new 
streamlined model helps facilitate our physical understanding of the instability, and 
also makes more tractable the analysis of its nonlinear dynamics. We concentrate on 
the linear theory in this paper, showing that the instability is restricted to a preferred 
range of intermediate wavenumbers and optical depths. We subsequently apply these 
general results to the inner B-ring and the C-ring and find that in both regions the 
ballistic transport instability should be near marginality, a fact that may have impor- 
tant consequences for its prevalence and nonlinear development. Owing to damping 
via self-gravity wakes, the instability should not be present in the A-ring. A following 
paper will explore the instability's nonlinear saturation and how it connects to the 
observed large-scale structure. 
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1 INTRODUCTION 

Like all solar-system bodies, the component particles of 
planetary rings must endure a hail of interplanetary me- 
teoroids impacting at speeds on the order of 10 km s~^. 
Collectively the most erosive projectiles lie in the size range 
of 10-^-10"^ cm, and it is estimated that Saturn sweeps up 
such particles at a rate > 10 kg s~^ (Durisen 1984; Ip 1984; 
Cuzzi & Durisen 1990, hereafter CD90). Hypervelocity im- 
pacts liberate significant amounts of material from ring par- 
ticles (some 10^-10^ times the mass of each impactor), and 
these ejecta re-accrete on to the ring, typically at a differ- 
ent radial location from where they started. The exchange 
of ejecta between nearby regions, referred to as 'ballistic 
transport', facilitates a redistribution of mass and angular 
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momentum on length-scales hh ~ 10 — 10 km and times 
~ 10^ - 10^ yrs that should control the large-scale evo- 
lution of Saturn's rings. Indeed, these scalings suggest that, 
over the age of the solar system, lO'^-lO* times the mass of 
the current rings has been transported in this way. 

Previous theoretical work shows that the ballistic trans- 
port process is adept at reshaping pre-existing gradients in 
surface density and composition. It can sharpen edges, such 
as those at the inner boundaries of Saturn's A and B-rings, 
generate the ramp features at the feet of those edges, and 
influence global colour gradients (Ip 1983; Lissauer 1984; 
Durisen 1984; Durisen et al. 1989, hereafter D89; Durisen et 
al. 1992; Cuzzi & Estrada 1998; Charnoz et al. 2009). But 
ballistic transport can also produce structure spontaneously 
from a homogeneous ring via a linear instability. This 'bal- 
listic transport instability' is thought to drive the 100-km 
wavetrains in the inner B-ring, and possibly the 100-km- 
wide plateaus and the low-amplitude 1000-km undulations 
in the C-ring (Durisen et al. 1992; Durisen 1995, hereafter 
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D95; Colwell et al. 2009; Charnoz et al. 2009). It is to the 
basic theory of baUistic transport and the hnear instabiUty 
that this paper will be devoted. 

The dynamics of ballistic transport has been success- 
fully described by a detailed global model constructed 
by Durisen and coworkers in the 1980s and 1990s (D89; 
CD90; Durisen et al. 1992; D95; Durisen et al. 1996; Cuzzi 
and Estrada 1998). The model incorporates a great many 
physical processes, but the mathematical formalism can 
be unwieldy and thus potentially obscure the fundamen- 
tal physics. In this paper we omit extraneous details and 
construct a simpler model that is easier to work with, yet 
remains sufficiently accurate. It is, in fact, almost identical 
to the leading-order Durisen formalism when expanded in 
the small parameter q = Vc/vc <^ I, where «o is the typical 
relative speed of the ejecta, and Vc is the orbital speed of 
the ring particles. (Note that Durisen and coworkers denote 
this ratio by x, which we reserve for a radial coordinate.) 
As a result, our model is local (the shearing sheet), and this 
makes the ejecta orbital dynamics easy to describe. The lo- 
cal formalism also permits the resulting transport terms of 
mass and angular momentum to be manipulated into one- 
dimensional integrals in convolution form. Being especially 
amenable to Fourier analysis, the linear theory is transpar- 
ent and nonlinear simulations are comparably straightfor- 
ward. Overall, the simpler formalism facilitates our physical 
understanding of the instability, and permits us to bring 
to bear the techniques of nonlinear dynamical systems. Ad- 
ditional physics can always be added later to sharpen the 
quantitative comparison with observations. 

In this paper we concentrate on the linear analysis of 
the system. It yields a simple instability criterion, which re- 
veals instability is facilitated, in particular, by the decrease 
in the ejecta emission rate near optical depths of 0.5. This 
drop corresponds to a transitional regime in which the disk 
becomes sufficiently dense that some of the liberated ejecta 
are reabsorbed by neighbouring particles rather than sent 
into orbit. In agreement with D95, we find that instability is 
suppressed at very low and high optical depths. But we also 
see that it is suppressed at very long and short wavelengths, 
and that unstable modes can propagate either radially in- 
wards or outwards, depending on wavenumber. Our analysis 
frames the problem in terms of two dimensionless param- 
eters: the mean optical depth tq and the 'ballistic Prandtl 
number' /i, which describes the relative efhciency of mass re- 
distribution caused by viscous stresses versus that caused by 
ballistic transport. Most of the uncertainties in the problem 
are packaged into n, which nevertheless can be tightly con- 
strained. In the A-ring /i is relatively large, because of strong 
self-gravity wakes, and therefore viscous diffusion smears out 
potentially unstable modes. In both the inner B-ring and 
the C-ring fj, is smaller and linear modes may grow. How- 
ever, in these two cases instability is near marginality, within 
the uncertainties, because ro is small or large, respectively. 
Though instability is still likely to occur, the fact that the 
system is near marginality will influence its nonlinear de- 
velopment, potentially leading to low-amplitude saturation 
or bistability. A weakly nonlinear analysis of these cases, 
together with fully nonlinear numerical simulations, will be 
presented in following work. 




Figure 1. Diagram showing the rationale and main features of 
the shearing sheet modeh A small patch of disk is isolated and 
treated as a Cartesian sheet subject to rotation and a shear flow. 
Terms arising from the cylindrical geometry are dropped and the 
coordinates x and y point in the radial and azimuthal direction 
respectively. 

The paper is organised as follows. In Section 2 we detail 
the mathematical formalism that we use to describe the bal- 
listic transport process, culminating in the governing equa- 
tion for optical depth, Eq. (22). There we also discuss the 
functional forms for the rate of emission, the probability of 
absorption and the ejecta throw distribution, which we draw 
from fits to the numerical calculations of CD90. The linear 
stability analysis follows in Section 3, in which we present 
growth rates and a general instability criterion with applica- 
tion to Saturn's A, B, and C-rings. We discuss these results 
and conclude in Section 4. 



2 MATHEMATICAL FORMALISM 

We consider a local model of a particulate ring, the shearing 
sheet (Fig. 1), which is a convenient representation of rela- 
tively small-scale dynamics in a differentially rotating disk 
(Goldreich & Lynden-Bell 1965). Instead of dealing with the 
entirety of a ring system, with its unconstrained global struc- 
ture and boundary conditions, we concentrate on a small 
patch of disk centred on a fiducial radius ro and moving 
on a circular orbit with angular velocity fio = ^{ro)- The 
small patch can then be described by a Cartesian coordinate 
system (a;, y, z), with x and y pointing in the radial and az- 
imuthal directions respectively, and its differential rotation 
represented by a combination of uniform rotation, floSz, 
and linear shear, v = —Soxey, where x = r — ro- The 
rate of orbital shear is So = —ro{dQ/dr)o and is equal to 
|f2() in a Keplerian disk. The shearing- sheet approximation 
introduces fractional errors of order A/ro, where A is the 
characteristic length-scale of the dynamics we want to de- 
scribe. Therefore, its descriptions of 100-km-long waves in 
the B-ring and plateaus in the C-ring yield an error ~ 10~^, 
whereas the very slow 1000-km undulations in the C-ring 
give 10^ . Similarly, the ballistic transport process can be 
described adequately by the shearing sheet, as the errors 
introduced scale as kh/ro ~ 10~^-10~^. 

Let a{x,t) be the surface mass density of the ring, as- 
sumed from the outset to be axisymmetric. We take the 
particle size distribution to be fixed and assume there exists 
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a unique relation between a and the normal optical thick- 
ness T. Mass conservation furnishes us with the following 
evolution equation for the ring: 

dtCT + dx{cTUx) =1 - J, (1) 

where I{x,t) is the rate at which mass is gained, per unit 
area, via ballistic transport from other radii, and J(x,t) is 
the rate at which mass is lost. The radial drift speed within 
the ring is denoted by Ux and is instigated by both viscous 
stresses and the ballistic transport of angular momentum. In 
the following subsections we derive a formalism that supplies 
us with convenient expressions for T, J and the local mass 
flux a Ux ■ 

2.1 Characterising the properties of ejecta 
emission and absorption 

2.1.1 Rate of emission, R 

The local rate at which mass is liberated from the ring per 
unit area as a result of meteoritic bombardment is the ero- 
sion rate R[a{x,t)]. This suggests a local gross erosion time 
te = (7/R[a], which corresponds to the time it would take 
for the ring to be completely destroyed, in the absence of 
ejecta recycling. At r = 1, an estimate for to is 

if) {mf^)y 

where the yield Y is the ratio of liberated mass in a moto- 
oroid impact to the mass of the impactor (D95). The value 
of Y depends sensitively on the strength and composition 
of the ring particle's surface, with 'softer' particles releas- 
ing more material and thus taking larger Y. Though the 
physical state of ring particles is poorly constrained, ter- 
restrial laboratory experiments give Y w 10^-10^ for hard 
water- ice/silicate targets struck by impactors travelling at 
speeds ~ 0.1-10 km s'^ (Lange & Ahrens 1987; Frisch 1992; 
Koschny & Griin 2001a). Note that these estimates per- 
tain only to cratering (non-disruptive) impacts, which are 
the most common outcome in Saturn's rings. Disruptive im- 
pacts produce higher effective yields, the cumulative effects 
of which may influence regions of the rings in which particles 
are smaller and hence more prone to destruction (Estrada 
and Durisen 2010). 

2.1.2 Throw distribution, f 

Liberated material travels in a slightly inclined and slightly 
eccentric Keplerian orbit that intersects the ring plane at the 

location of emission, x, and a distant location, x + ^. The 
radial distance ^ between these two locations depends on the 
ejection velocity, the statistics of which is summarised in a 
'throw distribution' /(^), normalised such that //($)d^ = 
1. (Unless otherwise specified, all integrals are carried out 
from — cx) to oo.) Thus /(^) d^ is the proportion of ejecta that 
travels a radial distance between ^ and ^+d^ during its orbit. 
The distribution should exhibit a characteristic length-scale, 
the 'throw length' Zth- By considering the dynamics of ejecta 
trajectories (Subsection 2.2.1; Appendix A), we set 

^--^-^-2x^o^(io^Ft)(io^)^'^^-' 



where g = Va/{roiT) and here Ve refers to the mean ejection 
speed. Thus Zth is the maximum throw possible for the given 
mean ejection speed Ve. The magnitude of Ve, like the yield 
Y, depends closely on the physical state of the ring parti- 
cle's surface, with 'softer' particles possibly taking 1 m s~^ 
and 'harder' particles taking 100 m s^^ (D89; Frisch 1992; 
Koschny & Griin 2001b). Consequently, the small parameter 
Q varies between 10~* and 10~^, while Zth lies between 10 
and 1000 km, which encompasses the scales of the phenom- 
ena we aim to describe. 

2.1.3 Probability of absorption, P 

During their orbit, ejecta may be reabsorbed by the ring at 
either the radius of omission x or at the distant intersection 
radius a: -|- ^. The probability that an ejectum is absorbed 
at a; -|- ^ we denote by P[(7{x + ^, t), a{x, t)]. The reason for 
allowing P to depend on the surface density at the emitting 
radius x, as well as that at the distant radius x + ^, is to 
account for extended excursions whereby the ejectum passes 
through the ring plane multiple times. This outcome is likely 
only in optically thin regions of the disk, and thus the second 
dependence may be dropped for optically thicker regions. 
Note that the ballistic transport of individual ejecta occurs 
on an orbital period ~ 10 hr, much less than the erosion 
time te. Hence the time required to execute multiple orbits 
can be neglected. 

2.2 Direct mass transport 

Now that we have introduced R, f and P we can construct 
the mass gain and loss terms I and J^. The rate at which 
mass is lost from a radius x is R[a{x, t)], while the proportion 
of this mass that travels to an annulus a distance ^ away 
and of thickness d^ is /(^) d^. Finally, the fraction of this 
mass that is ultimately absorbed by this annulus is P[a{x + 
^,t),a{x,t)]. Now if we sum over all such annuli we obtain 
the total loss rate 

J{x,t) = R[cr{x,t)] I P[a{x + ^,t),a{x,t)]f{OdC (4) 

A similar argument accounts for the total gain of mass at x 
from all neighbouring annuli: 

I{x,t) = j R\a{x-i,t)]P\a{x,t),a{x-i,t)]f{i)&t (5) 

The gain and loss integrals, though potentially compli- 
cated, benefit from being one-dimensional. The integrals in 
the formalism of D89, on the other hand, are three dimen- 
sional and cover the two emission angles of the ejecta as 
well as the speed of emission. These may be more accurate 
but arc awkward in analytic work, and costly in numerical 
simulations. Moreover, the three degrees of freedom in D89 
should map approximately to a single throw distance ^ in 
the case oi g <ti 1. In Appendix A we show in detail how the 
two formalisms join up. 

2.3 Angulcir momentum transport 

Both the ring's viscous stresses and the radial redistribu- 
tion of ejecta lead to angular momentum transport and a 
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consequent radial drift of material. In order to calculate the 
mass flux associated with the slow radial motion, we need 
to sketch out the trajectories of the ejecta in the shearing 
sheet. 



2.3.1 Ejecta trajectories in the shearing sheet 

The orbit of a representative ejectum [x{t),y{t),z{t)] in the 
shearing sheet obeys the Hill equations 

X - 2noy - ^VlIx = 0, (6) 
y + 20oi = 0, (7) 
+ n§2 = 0, (8) 

where an overdot indicates a time derivative. (Through- 
out this paper wc neglect non-Keplerian effects arising from 
planetary oblatcncss.) If the ejectum is thrown from the ring 
at t = from position (0, 0, 0) (without loss of generality) 
and with ejection velocity (w, it will undergo the fol- 

lowing trajectory: 

X = [« sin Qot -|- 2ti (1 — cos Qot)] , (9) 

y= ^[-2u{l-cosnot)+v{4smnot-30.ot)], (10) 
w 

z = — sinClot. (11) 
$2o 

The ejectum will first return to the ring plane half an orbit 
later, at t = tt/Oq, at position {4v,—4u — 3nv,0)/^lo and 
with velocity (— w, —7v, —w). If it is not reabsorbed at this 
point it will continue its orbit, meeting the ring plane again 
at t = 27r/r2o at position (0, — Gttu, 0)/no with its original 
velocity, {u, v,w). A representative trajectory is described in 
Fig. 2 in three dimensions and also projected in the {x, y) 
plane. The blue dots denote ring crossings at even integer 
multiples of t = 7r/f2, and red dots the ring crossings at odd 
multiples. 

In the shearing sheet context, the role of (specific) an- 
gular momentum is played by the (specific) canonical aa- 
imuthal momentum Py ~ y + 2^02;, a quantity that is con- 
served in the Hill equations, and which differs from the y- 
velocity because of the Coriolis force. For a particle in a cir- 
cular orbital motion at a fixed x we have y = —^Qox and 
thus Py — ^Qox. The uniform radial gradient of this quan- 
tity in the shearing sheet corresponds (apart from a factor 
of ro) to the local radial gradient of specific angular momen- 
tum among the family of circular Keplerian orbits. For the 
ejectum orbit considered above, Py = v, which can be seen 
from considering the conditions at t = at the moment the 
particle is launched. But this value corresponds to a circu- 
lar orbit located at x — 2v/flo, exactly halfway between its 
two ring-crossing radii. The radius of this orbit is described 
by the dashed black line in Fig. 2b. So when an ejectum is 
launched, it acquires a small amount (positive or negative) 
of angular momentum; subsequently, it oscillates about the 
radius associated with its new angular momentum and is 
absorbed at one or other radial extremum. Translating this 
result into the terminology of the previous subsection, we 
can then say that material emitted from x and with throw 
distance ^ has py = ^flo {x + C/2). 



2.3.2 Radial mass drifts due to ballistic transport and 
viscosity 

We now consider the evolution of the angular momentum 
analogue of the ring itself, crpy. The conservation law for 
this quantity in an axisymmetric ring is 

dt[(T{x,t)py]+da:Q=TC-C, (12) 

where AC and £ are the rate of gain and loss of angular 
momentum respectively at x due to ballistic transport, and 
Q is the local momentum flux density. It is equal to 

Q = (TUxPy —lixy, (13) 

where Ilxy is the xy component of the viscous stress tensor. 
As ring material follows circular orbits predominantly, we 
take Py = ^^lox from now. 

The two ballistic transport terms K. and £ are straight- 
forward to construct. Consider ejecta released at location x 
and absorbed at a; -|- ^ in an annulus of thickness d^. The 
(specific) angular momentum of such ejecta is ^ilo{x + $,/2). 
As before, the rate of its emission is R[a{x,t)], and the pro- 
portion that travels to the annulus is /(^) d^. Finally, the 
fraction that is absorbed at a; -|- ^ is P[a{x + ^,t),a{x,t)]. 
Summing over all the neighbouring annuli yields the rate of 
angular momentum loss at x: 

Cix,t) = R[a{x,t)] J P[a{x + ^,t),a{x,t)] f{0 

X ^ao{x + i/2)dC (14) 
The gain rate can be constructed in a similar way: 

TC{x, t) = j R[a{x - i, t)] P[a{x, t),a{x - ^, t)] /($) 

X ifio(a;-^/2)de (15) 

The nature of the viscous stress Uxy in cold and 
dense particulate rings, such as Saturn's, is nontrivial and 
comprises various components, each of which can deviate 
from the familiar Newtonian prescription (see e.g. Latter 
& Ogilvie 2006, 2008; Schmidt et al. 2009). In the in- 
ner rings viscous transport is dominated by the 'collisional 
stress' (Shukhman 1984; Araki & Tremaine 1986; Wisdom 
& Tremaine 1988), whereas stresses arising from self-gravity 
wakes prevail in the A-ring (Salo 1992; Daisaka et al. 2001, 
Yasui et al. 2012). The internal processes governing both op- 
erate on the orbital time-scale, which is much shorter than 
te. It hence makes sense to treat the stress in the diffusion 
approximation and to introduce an effective viscosity v. The 
viscous mass flux then becomes 

3 

Hx!/ = —-floi'a. (16) 

Generally, the viscosity i/ may be considered a function of 
cr or optical thickness r. For simplicity, we treat v here as a 
constant, although the linear theory presented in this paper 
is easily adapted, as we describe below, to the more realistic 
situation of a density-dependent viscosity. 

Finally, by subtracting from Eq. (12) ^Qox times the 
mass-conservation equation (1), we can obtain an expression 
for the mass flux density a that appears in the governing 
equation (1), 

aUx = -^{K: + £)-3d4ua). (17) 
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Figure 2. A representative ejectum trajectory begun at (0, 0, 0) and evolved forward for 4 ring-plane crossings according to Eqs (9)-(ll). 
The initial ejection speed is {u,v,w) = (0.1,0.1,0.1). The first panel describes the trajectory in three dimensions with the green grid 
indicating the ring plane. The second panel projects this trajectory on to the (x, y) plane. The dashed black line denotes the radius that 
corresponds to the ejectum's angular momentum. In both figures the axes have been stretched relative to each other in order to bring 
out more clearly the details of the motion. 



This expression introduces new transport integrals defined 
through 

4 

and 



(18) 



c^ — {c-\noxj). 

"0 



(19) 



Thus K, and C are identical to T and J but with an extra 
factor of ^ in their integrands. 

Our governing equation for the surface density is then 



dta=I-J + \d^{lC + £) + Zvdla. 



(20) 



As is well known, the viscous transport of angular momen- 
tum leads to a radial spreading of mass. The mass diffusion 
coefficient in a Keplerian ring resulting from a uniform kine- 
matic viscosity is 3/^. 



2.4 Governing dimensionless equation 

Once we specify R, P, f and v, we have all the ingredients 
to solve for the evolution of the ring. To simplify the fol- 
lowing calculations, Eq. (1) is non-dimensionalised. Time is 
scaled by tc and space by Zth. Surface density is scaled by the 
reference density ai. We then define the dynamical optical 
depth via 



r{x,t) = a{x,t)/ai, 



(21) 



and so ai is the density associated with t — 1. Note that 
the dynamical r can differ from the ring's physical or photo- 
metric optical depth measured by Cassini; this is especially 
the case when there exist self-gravity wakes (e.g. Salo and 
Karjalainan 2003, Porco et al. 2008, Robbins et al. 2010). 
That said, in the regimes relevant to the inner B-ring and 



the C-ring, the discrepancy is not severe and we treat the 
various optical depths as approximately equal. Finally, we 
scale R by ai/tc and / by l//th. Both R and P hereafter 
will be considered functions of r. 

The scaled evolution equation for r is 

dtT = T-J+\d.^{lC + C) + ^idlT, (22) 

where we have introduced the ratio of the mass diffusion due 
to viscosity to the 'ballistic diffusivity'^: 



(23) 



We emphasize that this 'ballistic diffusivity' is only a dimen- 
sional estimate, and will see below that ballistic transport 
does not in fact lead to a diffusion of mass in the conven- 
tional sense. The four ballistic transport integrals can be 
worked into the following forms: 



T{x,t)= J R[r{x-^,t)]P[r{x,t),r{x~^,t)]fiOd^, 

(24) 

J{x,t) = R[r{x,t)] J P[r{x + ^,t),r{x,t)]f{^)d^, (25) 

K:ix,t) = I R[T{x^^,t)]P[Tix,t),rix~U)]^f{0 dC, 

(26) 

C(x,t) = R[T(x,t)] J P[r{x + ^,t),r{x,t)]tf{OdC (27) 
If we make the additional approximation that P depends 



^ If I' depends on a we can generalise in the linear theory, 
by replacing u by d{i'<j)/d<T evaluated at the unperturbed surface 
density. 
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only on the optical depth at the distant, non-emitting, ra- 
dius, then we may suppress the second argument of this 
function. This also means that P can be taken outside the 
integrals in the expressions for X and IC. The four expressions 
can then be written in a compact way, in which * denotes a 
convolution integral with respect to x: 



1-- 
K. - 



--P-{R*g), 



J - 

C ■ 



■ R-{P- 



5). 



(28) 
(29) 



Here the function g is defined by g{£^) = £,/{£,), and the tilde 
denotes a reflection, so that /(^) — /(— ^ and g{^) — g{~£,)- 

2.4-1 The ratio of mass transport coefficients, fj, 

The only control parameter that appears in Eq. (22) is fi, 
which can adopt different values depending on the dominant 
mode of viscous transport. In the A-ring, we expect u to be 
monopolised by the action of self- gravity wakes. Analysis 
of density wave damping gives ~ 30 — 200 cm^ s"'^ (Tis- 
careno et al. 2007) , an estimate that agrees with direct mea- 
surements of f from A'^-body simulations of self-gravitating 
particles (Daisaka et al. 2001, Yasui et al. 2012). In the in- 
ner B-ring and the C-ring, however, we expect little or only 
moderate wake activity. In these cases, estimates from ki- 
netic theory and simulations give u ~ 0.1 cm'^ s~^ for the 
C-ring and v > 1.0 cm^ s~^ for the inner B-ring (Salo et 
al. 2001; Schmidt et al. 2009, Yasui et al. 2012). 

On the other hand, the ballistic diffusivity li-^/tc could 
be estimated by observing the ranges in its components from 
Eqs (2) and (3). If /th and tc are treated as independent, then 
the ballistic diffusivity is poorly constrained, its value rang- 
ing over four orders of magnitude: Ith/tc ~ 1 — 10* cm^ s~^. 

In fact, Ith and ie are correlated because of the approx- 
imate relation Vc oc y"^/^ (noted in Durisen et al. 1992). 
Such a relation would be consistent with a certain fraction 
of the impact energy being transferred to the ejecta. If we 
denote this fraction by e we obtain: 



2 

^ ^imp 



Yvl, (30) 

where «imp is a typical impact speed. Because Ith/'tc ccYv^, 
this leads to a tight bound on the ballistic diffusivity, sub- 
ject to an estimate of the transfer efficiency e. Previous nu- 
merical and experimental studies in a variety of materials, 
such as basalt, glass, gabbroic anorthosite, and powdery re- 
golith (but not ice), yield e ~ 0.1 in the relevant Vimp range 
(O'Keefe & Ahrens 1977; Hartmann 1985; Rashev & Ahrens 
2007) , and we adopt this as our fiducial value. We then find 



^1.3x W 

to 



lOkm/s 



105 km 



100 g/cm2 



cm /s, (31) 



and we expect it to take values of some 10 cm^ s~^ and 100 
cm^ s~^ in the C-ring and in the A and B-rings respectively. 
Substituting these estimates into (23) yields ~ 1 in the 
A-ring and ^ ~ 0.01 in the B and C-rings. 

The behaviour of the dimensionless system (22) is thus 
controlled by a single tightly constrained parameter, fi. But 
there still remains a broad spread in the physical length 




Durisen et al. (1989) 
Skew normal 
Off-centred Gaussian 




Figure 3. The first panel shows the normalised mass emission 
rate R (blue solid line), as approximated by D89 and given by 
Eq. (35), and the simple absorption probability p(i) (red dashed 
line), given by Eq. (32). The second panel shows our approxima- 
tion of the throw distribution / as a function of throw distance 
^ (black solid line) as calculated by D89 (see Appendix A). In 
addition we plot more convenient analytic estimates: the skew- 
normal distribution (red dashed line) with a = 1, = 1 ^ud 
/3 = —3.5; and the off-centred Gaussian (blue dot-dashed lines) 
with ^0 = 0.5 and a = 0.6. 



and time scales of the problem. Nonetheless Eqs (2) and 
(31) may replace (3), and therefore estimations of both the 
physical length Ith and timescales te of ballistic transport 
can be reduced to the important dependency on Y, the only 
poorly constrained parameter. For instance, setting Y — 10* 
and assuming inner B-ring densities gives ~ lO'' years and 
kh ~ 500 km. 
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2.4-2 Probability of absorption 

This and the next subsections present expressions for the 

functions P, f and R and discuss their form and the relevant 
physics in each case. We begin by describing the probabihty 
of absorption P. 

The fraction of ejecta that manages to pass through 
the ring plane on a given pass is exp(— Tsi,i), where Tsi,i is 
the slant optical thickness at the radius of interception. The 
slant optical thickness will always be equal or greater than 
the normal optical depth n, with the two related by Tsi,i = 
Ti/ cos 6inc where ^inc IS thc auglc of incidence of the cjccta in 
the frame of reference of thc receiving ring particles (D89). 
If wc next suppose that all thc ejecta that manage to pass 
through the ring at the distant radius are then reabsorbed 
at the original emitting radius we may assign 

P{n) = 1 - exp(-T,i,i). (32) 

Where necessary, we denote this expression by P*-^^ to dis- 
tinguish it from the more detailed version presented below. 

However, if the emitting region is of low optical depth, 
an appreciable amount of material may pass through and 
undergo another orbit and hence another opportunity to be 
accreted at the distant radius. The probability that ejecta 
is absorbed on the second try is 



[1 - exp(-Tsl,i)] exp(-Tsl,e - Tsl,i), 



(33) 



where Tsi_c is the slant optical depth of the ejecta at the emit- 
ting radius. We now sum over the formally infinite rmmber 
of circuits the ejecta may complete and obtain a convergent 
geometric series. We can then write down the total proba- 
bility of absorption: 



P(ri,Te 



1 - exp(-Tsi,i) 



1 — exp(— Tsi,, 



(34) 



Where necessary, we denote this more advanced model by 

p(2) 

To make further progress, wc need to account for the 
angles of incidence 6'inc that appear in formulas (32) and 
(34), which should vary depending on the geometry of the 
orbit. Following D95, we let cos Sine take a single 'average' 
value Tp. We thus set rsi,i — n/rp and Tbi^c ~ Tc/rp. As is 
argued in D95, Tp ~ 0.5 is a reasonably accurate approxima- 
tion for a typical set of trajectories, with the error worsening 
in lower optical depth regions. 



2.4-3 Rate of emission 

In CD90 a detailed formalism is constructed whereby the 
intensity and angular distribution of ejecta from a ring layer 
may be numerically calculated. The approach borrows much 
from the calculation of light-scattering from a layer of parti- 
cles, treating the diffuse incident intensity of meteoroids as 
similar to that of incident photons. Thc results of laboratory 
experiments are also used to infer the ejecta emission prop- 
erties, as functions of (a) the relative velocities of a spherical 
ring particle and an impacting meteoroid, and (b) the angle 
between the impact direction and the surface normal. The 
single particle scattering function may then be obtained by 
integrating over the ejecta contributions over all points on 



the spherical ring particle. The ejecta distribution function 
proceeds directly. Finally, though the meteoroid influx is as- 
sumed to be isotropic in the heliocentric reference frame, it 
is aberrated in the ring reference frame because of Saturn's 
and the ring particles' orbital motion. 

The results of these calculations give i? as a function 
of the optical depth at the emitting radius, and the ejecta 
distribution / as a function of the ejection velocity vector. 
The former may be approximated, following D95, by the 
dimensionless analytic form 



R{t) = 0.933 



1 + 



1 exp{-T/Ts) 



(35) 



where Tb = 0.28 is a parameter. This expression is normal- 
ized such that R{1) = 1. 

In Fig. 3a we plot i? as a function of r. At low r the 
emission rate is small because the impact rate is small. As r 
increases so do the number of impacts, and consequently R- 
The drop in 7? at ro ~ 0.5 we attribute to a transition from 
an optically thin regime, in which almost all ejecta from a 
given ring particle are thrown into orbit, to an optically thick 
regime, in which an increasing amount of liberated ejecta is 
reabsorbed immediately by neighbouring ring particles. At 
lower r it is possible for impact ejecta to leave the ring plane 
from both sides, whereas at higher r ejecta can leave from 
only one side because of the intervening particles. For suffi- 
ciently large r the emission rate relaxes to a constant value, 
as incoming meteoroids penetrate only to an optical depth 
of order unity and thus can only dislodge a fixed amount of 
material. The tendency for impact ejecta to be immediately 
accreted by neighbouring particles weakens R's dependence 
on r, relative to P's dependence (see Fig. 3a). This discrep- 
ancy plays an important role in the ring's stability and is 
especially marked near r = 0.5. 

2.4-4 Throw distribution 

Thc CD90 distribution function is not in the format required 
by our formalism as it depends on the ejection velocity: the 
two spherical angles that determine its orientation in the 
ring particle frame, 6 and 0, and its magnitude Vc- The dis- 
tribution function we have introduced, instead, depends on a 
single variable f , the throw distance. However, in the context 
of the transport integrals, an approximate mapping between 
{9, 4>, Vc) and ^ is possible, which allows us to translate the 
CD90 distribution to the / used in this paper. These de- 
tails are left to the Appendix. In Fig. 3b, we plot the / 
constructed in this way alongside two convenient analytic 
approximations, (a) the shifted skew-normal distribution 

fiO = exp[-(e - eo)V(2a')] (1 + erf (^ - 6)]) , 

(36) 

where erf is the error function, a is the 'standard deviation', 

P is thc 'shape parameter' which measures thc skewness, and 
the shift is ^o, and (b) an off-centred Gaussian, 



/(0 = ^^exp[-(C-Co)V(2a')] 



(37) 



with shift ^0 and standard deviation a. Each of these exhibits 
the main feature of the realistic distribution, which is its 
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characteristic asymmetry. Parameters that best match the 
computed / of D89 arc a = 1, = 1, and /3 = —3.5 for the 
skew normal and a = 0.6 and = 0.5 for the off-centred 
Gaussian. 

The asymmetry in / plays an important dynamical role, 
far greater than the asymmetry arising from the cylindrical 
effects that controlled the early studies of edge sharpening 
and redistribution (Ip 1983; Lissauer 1984; Durisen 1984). 
This outward bias is a consequence of meteoroids mainly 
striking the leading hemispheres of the ring particles. Ejecta 
are usually backscattered in these cratering impacts and 
thus adopt prograde orbits with v > 0, as in Fig. 2. Ring 
particles suffer more impacts on their leading faces because 
of aberration effects that follow from the orbital motion of 
the ring particles around Saturn and the motion of Saturn 
itself through the meteoroid flux. In particular, the orbital 
motion of the particles gives rise to a 'headwind' of matter 
that increases both the number of impacts on their leading 
faces and the impact velocities «imp. Terrestrial experiments 
on ice targets tell us that the ejecta yield obeys Y oc vf^^p 
(see CD90), and so the ejecta population (and hence its dis- 
tribution /) will be dominated by such prograde ejections. 



2.5 Integral relations 

Before applying our formalism to the question of the ring's 
stability, we sketch out some general results regarding the 
global conservation of mass, momentum, and energy. Note 
that the mass and angular momentum of impacting mete- 
orites have been neglected and thus do not appear in these 
balances. This is justified on the basis that Y ^ 1. The injec- 
tion of energy, however, may be more significant and could 
provide a minor source of ring particle velocity dispersion. 
Our formalism, however, does not treat the particles' ran- 
dom motion explicitly and this effect appears only through 
the kinematic viscosity f. See Durisen et al. (1996) for a 
more detailed discussion of this physics. 

Suppose Eq. (20) is solved for a ring of finite radial 
extent on an unbounded domain, meaning that either the 
density is of compact support, or decays sufficiently rapidly 
as X ^ ±oo. The following equations for the moments of the 
density distribution can be derived: 

^- J o-da; = 0, (38) 



dt . 

d_ 
di 



I' 



xdx = 0, 



f IT dx = f adx. 
dtj J 



(39) 



(40) 



These three equations can be understood as describing the 
evolution of the total mass, angular momentum and energy 
of the ring system. 

It is clear that Eq. (38) expresses the conservation of 
mass. To derive this equation, we note that 



Jdx = 0, 



(41) 



which is physically obvious and follows mathematically from 
a change of variables in one of the integrals. 



As we have already noted, Py = y + 2Clox plays the role 
of specific angular momentum in the local approximation. 
For a simple orbital motion in which x = 0, y = — ^ilox and 
z = 0, which is the local representation of a circular orbit in 
the reference plane, we have py — ^n^x. Therefore Eq. (39) 
expresses the conservation of angular momentum. To derive 
this equation, we note, using a similar change of variables, 
that 



jlxdx — jjxdx = j Kdx = j 



C dx. 



(42) 



and carry out an integration by parts for the viscous term. 
The specific energy in the local approximation is 



1/.2| .2, ■2\ 3 o2 2 I 1 o2 2 



(43) 



which equates to —^Qqx'^ for simple orbital motion. There- 
fore Eq. (40), which is derived in a similar way, states that 
the orbital energy of the ring decreases in time as a result of 
viscous dissipation. Ballistic transport, however, conserves 
mass, angular momentum and orbital energy. Eq. (40) also 
implies that the standard deviation of the mass distribution 
of a finite ring increases in time, i.e. that the ring spreads. 
This spreading is not (directly) counteracted by ballistic 
transport. Thus the edge-sharpening effects observed in sim- 
ulations should bo understood as a reshaping of a spreading 
ring feature and may be contrasted to situations where the 
spreading is actually halted, for instance by a shepherding 
satellite. In such cases Eq. (40) is modified by an external 
torque or angular momentum fiux reversal. 

The integral relations (39) and (40) do not generally 
hold when the equations are solved on a periodic domain 
because the angular momentum of material passing through 
the radial boundaries is not conserved. 



3 LINEAR STABILITY 
3.1 Dispersion relation 

Having set up a mathematical formalism to tackle the 
physics of ballistic transport, we now apply it to the question 
of a planetary ring's linear stability. We reserve the problems 
of ring edge sharpening and structure for future work. 

Consider a patch of disk in the homogeneous equilib- 
rium state of r = ro with fixed parameter. To case the 
exposition we assume for the moment that the absorption 
probability P depends only on the local optical depth at the 
emitting radius, i.e. P = P^^\ A disk of uniform r is an 
equilibrium because T = J and the local mass flux through 
the domain is uniform. The mass flux is in fact equal to 



PoJ?o ^ ^fd^ = -PoRo{0, 



(44) 



where a subscript indicates evaluation at To , and (^) is the 
expected (mean) throw distance. Note that, in the homoge- 
neous equilibrium we consider, there is no contribution to 
the radial drift velocity from viscous effects. Also note that 
if the throw distribution / is completely symmetric there is 
no drift velocity in the equilibrium state, as (^) = 0. Bal- 
listic transport then instigates no net angular momentum 
transport and hence no net radial drift. 
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We next superimpose on this state a small perturbation, 
f. After linearising Eq. (22), we obtain 

dtT = {RoPo - R'oPo)t + R'oPo{t * /) - RoPo{t * /) 
+ k^x [{RoPo + Pi>Po){Or + RoPo{t * g) + RoPoir * fl)] 
+ f^dlr, (45) 

where a prime indicates a derivative with respect to r. So- 
lutions exist of the form f oc e''*+'*'"^, where s is a (complex) 
growth rate and k a (real) wavenumber. When t has this 
form, the convolutions are easily evaluated as, e.g., 

f * / = F{k)f, (46) 

where 

F{k) = I fiOe-'^^dC (47) 

is the (non-unitary) Fourier transform of /(^). This result 
can be seen as a consequence of the convolution theorem, 
since the Fourier transform of t is proportional to 5{k). We 
thus obtain the dispersion relation 

s = RoPo - flo-Po + R^PoF{k) - RoPoF{-k) 

+ iife [{RoPo + RoPo){0 + R'oPoG{k) + RoPoG{-k)] 

- t^k\ (48) 

where G{k) = \F'{k) is the Fourier transform of g{£,), 
and the prime indicates differentiation with respect to the 
wavermmbcr k. Since /(^) is real, F{—k) = F{k), where the 
overline indicates complex conjugation. Noting further that 
F(0) = 1 and G(0) = {0, we find 

s = RqPo H{k) - RoPi, Wk) - MA;^ (49) 

where 

H{k) = F{k) -I- Ik [F'{k) + F'(0)] . (50) 

The first two terms in Eq. (50) come from I and J (ballistic 
mass transport) whereas the bracketed terms come from AC 
and £■ (ballistic angular momentum transport). Note that 
the real part of H comes from the even part of /, while the 
imaginary part of H comes from the odd part of /. Thus 
linear modes will manifest as travelling waves unless the 
distribution / is completely symmetric, as noted in D95. 

3.2 Stability criterion 

As a consequence of Eq. (49), the real part of s is simply 

Re(s) = {R'oPo - Ra)P») Re [H{k)] - fik\ (51) 

This isolates in a neat mathematical way the various in- 
gredients governing the mode's potential growth. The first 
bracketed factor in the first term accounts for the combined 
effect of the emission and absorption rates of the ejecta; the 
second factor summarises the influence of the throw distri- 
bution; and the last term introduces viscous damping. It is 
the first term that is responsible for instability, and it must 
arise from either the form of the absorption/emission profile, 
the peculiarities of the throw distribution, or a combination 
of the two. 

'Realistic' distribution functions /, such as those given 



in Section 2.3.3 and Fig. 3b, yield Re[ii"(fc)] < for all k. In 
fact, the real part of H can be positive only for exceedingly 
asymmetric and/or skewed distributions such as the delta 
function (single throw distance) considered by D95. It fol- 
lows that instability is controlled by the first factor in the 
first term in Eq. (51). As a consequence, we can immediately 
derive a necessary condition for instability: 

dlnP dlnR 

This condition is satisfied when the dashed curve in Fig. 3a 

is steeper than the solid curve (since the plot is logarithmic). 
So in order to obtain a growing mode, the rate of change of 
absorption must outstrip the rate of emission as we increase 
T. This mates intuitive sense. Consider a small overdensity 
upon a uniform ring: as a result of the local increase in r, 
both the absorption P and emission R of ejecta will adjust. 
If inequality (52) holds then the ring will absorb more ejecta 
in relative terms than it releases. As a consequence, material 
will start building up at that point and the overdensity will 
increase, which in turn will aid the accumulation of even 
more mass, and so on. Similarly, in an underdensc portion 
of an otherwise uniform ring, emission will out-compete ab- 
sorption and the underdensity will be exacerbated. The drop 
in R near r = 0.5, witnessed in Fig. 3a, almost guarantees 
instability for t > 0.5 because dP/dr > 0. 

The stability criterion (52) is only a necessary condi- 
tion because we have yet to include the damping effect of 
viscosity. A sufficient criterion for instability, involving both 
/i and r, is presented in Section 3.3.2. We can, however, 
establish some results in the limit of small and large k. At 
very small scales, i.e. large k, the Riemann-Lebesgue lemma 
tells us that both F{k) and F'{k) go to zero for realistic /; 
they do so exponentially fast if / is infinitely differentiable. 
Hence the —fik^ term in Eq. (51) dominates in this limit and 
modes are viscously damped, as is expected. On the other 
hand, on long wavelengths for which k is small, a Taylor 
series expansion reveals that Re[H(k)] oc fe*. Therefore vis- 
cous damping (oc —k^) will dominate in this limit as well, 
and very long modes will decay. This rcfiects the fact that 
ballistic transport is weak far beyond its throw length. In 
summary, if the ring is unstable, growing modes are limited 
to a band of intermediate wavelengths, as in convection or 
gravitational instability. In the next subsections these results 
will be illustrated numerically. 

In fact the Taylor expansion of H{k) is 

H{k) = ~F"'{0)k' - ^F""{0)k' + 0{k% (53) 

in which F"'(0) = i(C^) is imaginary and F""{0) = (?'') 
is real and positive. Terms proportional to k or k^ are ab- 
sent because of cancellation between ballistic mass transport 
and ballistic angular momentum transport; this result is re- 
lated to the fact that ballistic transport conserves the first 
and second moments of the surface density distribution (cf. 
Section 2.4). The absence of a k^ term shows that there is 
no 'ballistic diffusion' as such, although there are dispersion 
and hyperdiffusion. 

Before continuing on to the direct computation of 
growth rates, we examine the stability when the more ad- 
vanced probability absorption model in Subsection 2.3.1 is 
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adopted, and P — P^^\ Now P depends on emitting optical 
depth Tc, as well as the optical depth at the intercepting ra- 
dius Ti: P = P''^' (ti, Te). The linear analysis can be worked 
through as before, and we find the real part of the growth 
rate is: 



Re(s) = 



R'qPo — Ra ( — 



+ Ro { jT- 



Re [H{k)]-fik'^ 



(54) 

Here a subscript indicates evaluation at ri = Tc = tq. The 
necessary condition for instability, analogous to (52) , is then 



ainP dlnR dlnP 
> — ; h 



OTi ar OTc 

Because P is a decreasing function of Tc, from Eq. (34), the 
additional term on the right-hand side will be negative. As a 
result, the instability criterion will be easier to satisfy than 
(52). Note that the term in large square brackets in Eq. (54) 
is proportional to —{A — B) in D95. 

3.3 Specific examples 

3.3.1 Growth rates 

We now employ the choices of P, R and / introduced in 
Section 2.3 and compute growth rates explicitly. The distri- 
bution function we set equal to either the off-centred Gaus- 
sian or the skew-normal distribution, and we mainly take 
P — P'^' for simplicity. Unfortunately, the Fourier trans- 
form of the skew normal (36) is an intractable convolution. 
But the transform of the off-centred Gaussian (37) takes a 
simple form. In this case 



F{k) = exp (- 
and the real part of H{k) is 
ReiH) = -1-f [(1 



1 2,2 
-2« ^ 



ifc^o) 



^a^fc^) cos(fc^o) 



+ ^k^o sin(fc5o)] e 



''fc^/2 



(56) 



(57) 



This is negative for all fc 7^ provided that S,o/a. < 4.1595. In 
other words, Re(//) < unless the distribution / is narrowly 
confined to a single throw length ^0. The best fit to the 
D89 distribution gives ^0/1 ~ 0.83, which is well within this 
limit. The single throw length distribution of ^/a 2> 1 was 
examined in D95. In this case H is positive and hence the 
distribution can drive instability independently of the P and 
R profiles. As a result, the dispersion relation that ensues is 
more complicated (but unrealistic). 

In Fig. 4, we plot the real part of the growth rate ver- 
sus the wavenumber k when tq — 0.5 for two values of /i 
and for two different distribution functions /. The skew- 
normal distribution is represented by the red dashed curve 
and the off-centred Gaussian by the blue solid curve. For a 
lower value of /i, equal to 0.03, both curves exhibit growth 
on a characteristic range of wavenumber. Roughly, there is 
no growth on scales less than a throw length Ith, nor on 
scales longer than about 6Zth. The fastest growth occurs 
on a scale near (3/2)Zth with e-folding time ~ Ate. As ex- 
plained earlier, viscous diffusion dominates both long and 
short scales, because ballistic transport is inefficient in each 
limit. Note, however, that the skew-normal growth curve ex- 
hibits growth on slightly shorter lengths. Otherwise there is 
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Figure 4. Real growth rates versus wavenumber k for two values 
of pu and for two models of the distribution function /. The solid 
blue curve represents the off-centred Gaussian of (37) with go = 
0.5 and a = 0.6. The dashed red curve represents the skew-normal 
distribution (36) with a = go = 1 and fi = —3. In all cases 
the background optical depth is r = 0.5. The more viscous case 
fi = 0.15 does not exhibit instability at any k, but for fi = 0.03 
instability emerges on an interval of intermediate wavenumber. 




Figure 5. Real growth rates versus k for various values of tq. The 
off-centred Gaussian is used (as in in Fig. 4). The fi parameter 
is held fixed and is equal to 0.006. Note that instability occurs 
only on an intermediate range of optical depth t; both small- and 
large-T rings are stable. 



good qualitative agreement between the two distributions. 
For larger /x, growth is extinguished as viscous diffusion be- 
comes more efficient than ballistic transport and the wave 
modes are smoothed out before they can grow. Note that the 
growth rates (and their behaviour) are quantitatively con- 
sistent with the growth rates computed in D95 (his Fig. 8). 
In Fig. 5, growth rates are given for different tq at a 
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fixed fi. Here only the ofF-centred Gaussian lias been em- 
ployed. The figure illustrates the fact that growth occurs 
only on intermediate optical depths, around tq < 1: very 
large and very low optical depths do not exhibit instability. 
This reflects the relative profiles of the emission and ab- 
sorption functions R and P given in Eqs (32) and (35). In 
particular, instability favours the optical depths where -R's 
rate of change with r decreases and becomes negative. 

3.3.2 Stability criterion 

The stability criteria (52) and (55) are only necessary condi- 
tions because viscous diffusion has been omitted. A sufficient 
condition for instability must incorporate its influence and 
thus involve the parameter ^ as well as tq. From the form 
of the dispersion relation, marginal stability occurs when 
Re(s) = dRe(s)/dfc = 0. Given tq, we thus have two equa- 
tions to solve for a unique critical wavenumber k and critical 
^. This is accomplished numerically, and we plot the ensu- 
ing marginal stability curve /i = fJ.{To) in the 2D parameter 
space of (/i, To). 

In Fig. 6 we present the marginal stability curves for 
the two models of the absorption function P, given by (32) 
and (34). The distribution function / takes the off-centred 
Gaussian profile. Regions below the curves are unstable and 
regions above the curves are stable. The red dashed curve 
represents the more realistic P^'^\ which incorporates the 
variation in r at both the emission and intersection radii. 
The blue solid curve represents the simpler P^^\ which ac- 
counts only for r at the intersection radius. As is clear, at 
moderate to large r > 1 the two curves are much the same; 
this is because most ejecta in more optically thick regions 
undergo only one half or one orbit before being reabsorbed. 
The r at the emitting radius is then unimportant. However, 
at low r the two curves deviate, because ejecta undergo mul- 
tiple orbits more frequently, an outcome that is not modelled 
by the simpler P'^' model. The net effect of these multiple 
orbits is to enhance instability. So for a given low tq the 
critical /i can be double that predicted by the simple model. 

Perhaps the most important result here is that the max- 
imum /i that permits instability is remarkably low (~ 0.08). 
Viscous diffusion must be much less than ballistic transport 
or else the instability is washed away. Given our estimates 
for n, such a low value immediately rules out the A-ring 
as a venue for the ballistic transport instability: we find /i 
is ~ 1 in the A-ring. On the other hand, the estimates on 
H for the B and C-ring (~ 0.01) suggest that instability is 
possible, but only barely. From Fig. 6, the critical /i is 0.032 
when To = 1.1 and 0.042 when t = 0.1. If instability occurs 
in these regions then it may be near criticality, a fact that 
should influence its nonlinear saturation in important ways. 
We discuss this issue further in Section 4. Finally, we note 
that in the outer B-ring, where tq > 3, the critical jj, is tiny, 
and instability suppressed. 

Generally, for a given (sufficiently small) /i there exists 
an interval of r in which instability occurs (cf. Fig. 5, and 
D95). Instability is suppressed at both high and low r. At 
high optical thicknesses, the mechanism of instability be- 
comes weak, as both P and R have similar dependences on 
r (cf. Fig. 3a). Over and under-densities are only mildly ex- 
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^0 

Figure 6. Curves of marginal stability in the (tq,^) plane. The 
region below the curve is unstable and the region above the curve 
is stable. The solid blue curve corresponds to the case where the 
simple 'one orbit' model P^^^ is employed, Eq. (32). The dashed 
curve represents the case where the more realistic 'multi-orbit' 
model P(-^' is employed, Eq. (34). 



acerbated, and potentially unstable modes grow too slowly 
to escape viscous diffusion. For similar reasons, the instabil- 
ity mechanism weakens at very low r. In this limit as well, 
both P and R vary similarly with r and the first term in 
(51) is small as a result. Viscous damping again overwhelms 
potentially growing modes. 

These results are consistent with the stability bounds 
computed in D95. If the critical effective yields Y in D95 
are translated to critical then the curves in Fig. 7 agree 
to within a factor of 2 (R. Durisen, private communication). 
This is encouraging agreement given the different distribu- 
tion functions used (see the Appendix). 

3.3.3 Phase speeds 

Generally speaking, unstable modes manifest as travelling 
waves, because s is complex. Only if the distribution / is 
perfectly symmetric will modes grow in place. We define a 
mode's phase speed Cp via 

Cp = -ilm(s) = - i {RPy, Im[i/(fc)], (58) 

where we have used the simple model for P. In Fig. 7 we 
plot both the phase speed (red solid curve) and the real part 
of the growth rate (blue dashed curve) versus wavenumber 
k for the off-centred Gaussian model of /. 

A striking feature of Cp(fc) is that it changes sign as k 
increases. Longer growing modes propagate outwards, while 
shorter growing modes propagate inwards, leaving a critical 
k at which a mode grows monotonically. Thus the direction 
of propagation is not strictly tied to the bias in the throw 
distribution /. Having said that, most modes (including the 
fastest growing) travel radially inward. In addition, as k be- 
comes large, Cp asymptotes to a constant value. 
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Figure 7. Phase speed Cp and real growth rate as functions of fc, 
for r = 0.5, = 0.01. The off-centred Gaussian is employed for 



The sign reversal in the phase speed can be observed 
mathematically from the following argument. In the limit 
A: — >■ 0, we have H{k) ~ —j^i{(,^)k^ and so Cp ~ 
j^(-RP)o{^^)fe^. In the limit fc oo, we have H{k) ~ 
-iF'(0)fe = \i{i)k and so Cp \{RP)'o{Q. The rele- 
vant throw distributions are skewed such that both (5) and 
(5^) are positive, so a sign reversal in Cp must occur at in- 
termediate k. 

This behaviour can also be understood in more physical 
terms. Consider very short modes fc 3> 1 with wavelengths 
much less than both /th and the throw length's standard 
deviation (a for the off-centred Gaussian). In such a limit, 
the influence of the many small undulations are 'averaged 
away'. As a result, the mass transport imparts nothing to 
the collective motion. However, there will be a net angular 
momentum transport which will excite a net radial flux of 
material. With an outward bias to / this leads to an inward 
drift. 

Consider now a mode with wavelength much longer 
than Zth- Ballistic transport is then limited to relatively short 
distances. Suppose that mass is thrown almost entirely out- 
ward and that higher-density regions emit more mass. A 
density minimum will increase because it is receiving more 
material from the denser disk inwards to it than it can emit. 
Conversely, a density maximum will decrease, because it is 
emitting more mass than it is receiving from the (less dense) 
disk inwards to it. As a consequence of this differential mass 
transport, the entire wave-form will crawl outward, in the 
same direction as the throw asymmetry. 

The propagation speed of the fastest growing mode is 
~ 0.2 Ith/te which lies between 0.1 and 10 m yr~^. The waves 
move particularly slowly. This is consistent with Voyager 
and Cassini observations of inner B-ring structure that sug- 
gest these undulations have not travelled appreciably over a 
30 year period (Colwell et al. 2009). According to the above 
estimate, a linear mode will have propagated between a few 
metres and a few hundreds of metres in that time. The upper 



limit is just within the range of detection, but the lower limit 
certainly is not. The general consistency is encouraging and 
suggests a potential constraint on the typical ballistic trans- 
port speed hh/te- It must be stressed, however, that the 
observed structures are most likely nonlinear waves which 
may propagate at a different speed from the linear modes. 



4 CONCLUSION 

In this paper we construct a mathematical formalism that 
describes the ballistic transport in planetary rings. Tak- 
ing advantage of the relative smallness of the characteristic 
throw length Ith, which is connected directly to the param- 
eter Q <^ 1, we employ the local shearing sheet model and 
manipulate the transport terms into simple ID integrals in 
convolution form. The resulting main equation is simple to 
work with, both analytically and numerically, and its results 
easy to interpret. Moreover, it is almost as accurate as the 
classic formalism of D89, with relative errors probably of 
order g ~ 10"^ - 10"^ In this paper we deploy the model 
on one facet of the ballistic transport process: the linear 
instability (D95). But it may also be used to study ramp 
formation and steep edges. 

We derive the linear theory of the instability and ap- 
ply our results to Saturn's A, B, and C-rings. The stability 
analysis can be framed conveniently in terms of two param- 
eters: To the background optical depth and /i the ratio of 
mass transport coefficients due to viscous diffusion and bal- 
listic transport. We find that, for realistic profiles of the 
absorption probability P, the ejecta emission rate R, and 
the ejecta throw distribution / (D89, CD90), instability is 
pervasive for low and intermediate r. Actually, instability 
relies on the fact that R increases more weakly with r than 
P does, and it is especially exacerbated by the drop in R for 
r > 0.5 (cf. Fig. 3a). Near this optical depth, a small increase 
in density instigates a fall in emission, and a concurrent rise 
in absorption. The small overdensity is hence reinforced and 
the process runs away. The drop in R at this to we attribute 
to a transition from an optically thin regime in which almost 
all ejecta from a given ring particle are thrown into orbit, 
to an optically thick regime in which increasing amount of 
liberated ejecta is reabsorbed immediately by neighbouring 
ring particles. 

However, the ballistic transport instability is vulnerable 
to viscous diffusion which can smear out potentially grow- 
ing modes. In fact, the critical ^ above which instability is 
always extinguished is low: fic ~ 0.08. In the A-ring, self- 
gravity wakes dominate and viscous diffusion is relatively 
efficient; here fj, takes values ~ 1 and instability never oc- 
curs. In the B-ring, n ~ 0.01, which precludes instability 
in its extremely dense outer regions, where r > 3. For these 
large optical depths, Fig. 6 indicates fic is tiny^. On the other 
hand, the inner regions of the B-ring exhibit ro ~ 1 with a 
corresponding critical value of /Xc ~ 0.032. Thus instability 

^ The large-scale 100-km structure observed in the outer ring is 
most Ukely generated by something other than ballistic transport, 
perhaps electromagnetic instability (Goertz and Morfill 1988, 
Shan and Goertz 1991). 
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can occur but should be near marginality. The case is the 
same in the C-ring, where at tq — 0.1 we have /Xc ~ 0.042. 

While we can tightly constrain the governing stability 
parameters, jj, and mean optical depth to, the physical length 
and time scales of the phenomena arc less easy to tie down, 
owing to uncertainties in the physical state of ring parti- 
cle surfaces (and hence Y). Thus an exact determination of 
the preferred length-scales of the structures cannot yet be 
made. However, we may be able to better constrain these 
by matching the results of nonlinear simulations with the 
observed structures in detail. 

The fact that instability may be near marginality in 
both the B and C-rings may have important dynamical con- 
sequences. On the one hand, instability may saturate at a 
low amplitude. Indeed the C-ring exhibits gentle undulations 
of 1000-km wavelength. The C-ring, however, also supports 
100-km plateau structures of relatively large amplitudes at 
slightly larger radii (Colwell et al. 2009). Can the ballistic 
transport instability generate both sets of structures con- 
currently? This might suggest that ring properties in the 
C-ring vary rapidly with radius, yielding Zth ~ 1000 km 
near r = 80000 km and Zth f« 100 km at r = 90000 km. 
But it is also possible that the nonlinear dynamics permits 
the coexistence of both dominant lengthscales; if so, detailed 
numerical simulations may help in showing how. 

On the other hand, a system near marginality could ex- 
hibit bistability, whereby both a linearly stable homogeneous 
state and an active large-amplitude state are supported. In 
fact, the inner B-ring displays adjoining flat zones and active 
'wave zones' between r = 93000 km and r = 98000 km (Col- 
well et al. 2009). Can we associate these two regions with 
the inactive (homogeneous) and active (wavelike) states of 
a bistable system? If so, what determines the arrangement 
of the states; do the states propagate into each other at a 
'front'; and at what speed would this occur? Alternatively, 
the ring viscosity may vary sufficiently over the inner B- 
ring to permit instability in some regions and not in oth- 
ers. The main questions are then: what is the cause of the 
large-scale viscosity variation and how much would it need 
to change to explain the observations? Because the system 
is near marginality here, the viscosity need not change by a 
great amount. 

Some of these questions we hope to answer in future 
work where we will present the nonlinear theory and simula- 
tions of the instability. There we will exhibit the formation 
of wavetrain solutions, their stability to secondary modes, 
and the dynamics of potentially bistable regions of the disk. 
These various behaviours and patterns we will subsequently 
connect to observations. 
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APPENDIX A: CONNECTION WITH THE D89 
FORMALISM 

Most of the theoretical apparatus of D89 can be massaged 
into the form we present here by taking the local approxi- 
mation, expanding in the small parameter g and retaining 
only leading-order terms. A number of additional assump- 
tions then need to be made in order to connect the loss and 
gain integrals of D89 with the simpler versions that appear 
in Section 2. We sketch out these details now. 

In the formalism of D89 the ballistic transport terms 
are generally of the form 

poo r2-K pi 

A= / / / e{m^i)g{Q,(j),e)dQd(t>dcose, (Al) 
Jo Jo J -1 

where the angles (f> and 9 delineate the ejecta velocity ori- 
entation and g its magnitude (in units of rQo)- The angle 
6 denotes the angle between the y axis and the ejecta ve- 
locity, and (i) the angle between z axis and the projection of 
the velocity on to the {x, z) plane (see Fig. 2 in D89). The 
function is either P or R, and g is the CD90 ejection dis- 
tribution function. Recall that Tsia is the slant optical depth 
of the ring at the distant radius of ejecta re-intersection. The 
potential dependence of P on rsi.e is not denoted for ease of 
presentation. 

First we make the approximation rsi,i = ri/rp, that is 
we treat the slant optical depth as a simple function of nor- 
mal optical depth only, as assumed in Section 2.3.1. We then 
can set n = T(ri), where n is the radius of the distant rein- 
tersection point. This radius may be expressed in terms of 
the ejecting radius Te and (g, (/>, 0): 

_ [" 1 -h 2g cos g -h (cos^ e + SVC? e cos^ 0) 1 , . , 

^' ~ [ 1 - 2^. cos 61 - (cos2 e + sir? 6 cos^ 0) J ^ ' 

(sec D89 and Durisen et al. 1996). Expanding in small g 

gives simply 

n = re {1 + 4g cos e + 0{g^)) . (A3) 

Thus if we retain only terms up to g then n = n{g,6) and 
hence Q — Q{g,9) and the <j} dependence vanishes. All that 
matters now is the component of the ejecta velocity in the 
2/-direction. 

Next we define the throw distance 

i, = ri — rc — 4:grc cos 6, 

transform the 9 integral into a ^ integral, and swap the order 
of integration. This yields the form familiar from Section 2, 

/oo 
e[ri(ri+0]/(C)dC, (A4) 
- oc 

where we have introduced the new distribution function 

m)=rr ^d<t>dQ. (A5) 

Jo J|«|/(4r-e) ^STe 

D89 introduce an analytic approximation to the full nu- 
merical solution of g. It takes the form 

g = \sin(l)\^'*9exp{-9^-''^)h{g), (A6) 

where /i is a yet to be specified function of g. Using this 
estimate, the (f) integral in (A5) can be done immediately, 
leaving the sole g integral. The simplest choice for /i is a 



delta function 5{g — gc), for fixed gc- This means that all 
ejecta are released with the sauio speed. This then accounts 
for the g integral and we obtain 

/ cx cos-'(e//th)exp[-(cos-'(CAth))'-'], (A7) 

where we have set ?th = '^gcVe- A more realistic choice but 
with the same basic properties is 

h{g) = geM-{Q - Qcf/{2al)\, (A8) 

for some (narrow) dispersion . The g integral in ( A5) then 
must be accomplished numerically. The resulting function is 
plotted in Fig. 3b with od = 0.2. Note that in D95, g is 
mainly taken to be a distribution uniform over hemispheres 
with h a truncated power law, derived from hypervelocity 
experimental data (Durisen et al. 1992). 
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